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Abstract 

The paper presents the solutions for the zonal electrophoresis equations are obtained by analytical 
and numerical methods. The method proposed by the authors is used. This method allows to 
reduce the Cauchy problem for two hyperbolic quasilinear PDE’s to the Cauchy problem for ODE’s. 
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equations. The method is effectively applicable in all cases when the explicit expression for the 
Riemann-Green function of some linear second order PDE, resulting from the use of the hodograph 
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of constructing a multi-valued solutions. Compared with the previous authors paper, in which, 
in particular, the shallow water equations are studied, here we investigate the case when the 
Riemann-Green function can be represent as the sum of the terms each of them is a product of 
two multipliers depended on different variables. The numerical results for zonal electrophoresis 
equations are presented. For computing the different initial data (periodic, wave packet, the 
Gaussian distribution) are used. 
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I. INTRODUCTION 


In previous paper [1] , the efficient numerical method, allowing to get solutions, including 
multivalued 1 , are proposed in the case of the Cauchy problem for two hyperbolic quasilinear 
PDE’s. This method is based on the results of the paper [2] in which the hodograph method 
based on conservation laws for the two quasilinear hyperbolic PDE’s is presented. For the 
determination of the densities and fluxes of some conservation laws the linear hyperbolic 
second order PDE is used. As shown in [2], the solution of the original equations can easily 
be written in implicit analytical form if there is an analytical expression for the Riemann- 
Green function of mentioned linear hyperbolic equation. 

The paper [1] shows that one can not only write the solution in implicit analytical form, 
but also construct the efficient numerical method of the Cauchy problem integration. Using 
minor modifications of the results of [2] it is able to reduce the Cauchy problem for two 
quasilinear PDE’s to the Cauchy problem for ODE’s. From the authors point of view, the 
solution of the Cauchy problem for ODE’s, in particular, numerical solution, is much easier 
than the solution of nonlinear transcendental equations that must be solved when there is 
an implicit solution of the original problem. 

A key role for the proposed method plays the possibility of constructing an explicit ex¬ 
pression for the Riemann-Green function of the corresponding linear equation. This, of 
course, limits the application of the method. In fact, the number of equations to which the 
method is applicable is large enough. These include the shallow water equations (see, for 
example, [3, 4]), the gas dynamics equations for a polytropic gas [3, 4], the soliton gas equa¬ 
tions [4, 5] (or Born-Infeld equation), the chromatography equations for classical isotherms 
[3, 6, 7], the isotachophoresis and the zonal electrophoresis equations [8-12]. A large number 
of equations, for which the explicit expression for the Riemann-Green functions are known, 
are presented, in particular, in [2]. Classification of equations that allow explicit expressions 
for the Riemann-Green functions, is contained in [13-15] (see also [16, 17]). 

This paper presents analytical and numerical solution of the Cauchy problem for the 
zonal electrophoresis equations [8-12]. The choice of these problem, in particular, due to 
the fact that there is the Riemann-Green function which can be represented as a finite sum 

1 In [1] the solutions of the shallow water equations describing the braking waves are presented. 
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of two multipliers each of them depends on the different variable 


n 


$(i?\ R 2 |r\ r 2 ) = V k {R\ R 2 )Q k {r\r 2 ). 


( 1 . 1 ) 


k= 1 


It is shown below, that this type of function allows to significantly simplify the construction 
of solutions. 

Pay attention to the fact that in some sense, the proposed method is ‘exact’. Its real¬ 
ization does not require any approximation of the original hyperbolic PDE’s, which uses 
of the finite-difference methods, finite element method, finite volume method, the Riemann 
solver, etc. Also there is no need to introduce an artificial viscosity 2 . In other words, the 
original problem is solved without any approximations or modifications. The accuracy of 
the solution is determined by only the accuracy of the ODE’s numerical solution method. 

The paper is organized as follows. In Secs. II we repeat the slightly modified (and 
simplihed) results of the paper [2]. In Sec. II we formulate the problem for the zonal 
electrophoresis equations. Also we construct the Cauchy ODE’s problem which allows to 
obtain the solution of the original problem on the isochrones. In this section we present the 
results of calculating for the different initial data. 

II. REDUCTION OF THE CAUCHY PROBLEM FOR TWO HYPERBOLIC 
QUASILINEAR PDE’S TO THE CAUCHY PROBLEM FOR ODE’S 

Referring for details to [1, 2], here we give only a brief description of the method which al¬ 
lows to reduce the Cauchy problem for two hyperbolic quasilinear PDE’s to Cauchy problem 
for ODE’s. 

A. The Riemann invariants 

Let for a system of two hyperbolic PDE’s, written in the Riemann invariants R 1 (x,t), 
R 2 (x,t), we have the Cauchy problem at t — t 0 


R 1 + A \R\ R 2 )R l x = 0, R 2 + A 2 (R 1 , R 2 )R 2 x = 0, 


( 2 . 1 ) 


R l (j, t 0 ) = flj(i), R 2 (x, t 0 ) = -RoM. 


( 2 . 2 ) 


2 The effect of the grid viscosity does not occur due to the absence of approximation 
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where Rq(x), Rq(x) are the functions determined on some interval of the axis x (possibly 
infinite), A 1 (i? 1 ,i? 2 ), A 2 (R 1 ,R 2 ) are the given characteristic directions. 


B. Hodograph method 

Using the hodograph method for some conservation law (pt J r4’ x = 0, where <-p{R}, R 2 ) is 
the density, ^(i? 1 ,!? 2 ) is the flux, we write the equation [2] 


^r 1 # 2 + A(R , -R )$ri + B(R , R )$ij2 — 0, 
A(R\ R 2 ) = B(R\ R 2 ) = X2r1 


A 1 - A 2 ' 


A 1 - A 2 ' 


(2.3) 

(2.4) 


C. The Riemann—Green function 


Let the function $(R X , R 2 \r l , r 2 ) be the Riemann- Green function for equation (2.3). The 
function $(R X , R^r 1 , r 2 ) of variables R 1 , R 2 satisfies the given equation, and the function 
$(i? 1 ,i? 2 |r 1 ,r 2 ) of variables r 1 , r 2 is the solution of the conjugate problem 


$,. 1,2 - (A{r\ r 2 )$) r i - (^(r 1 , r 2 )$) r 2 = 0, 

($, 2 - A$) | rl=fll = 0, ($,i - 5 $)\ r2=R2 = 0, 


<T> =1 

I r 1 =R 1 ,r 2 =R 2 


(2.5) 

( 2 . 6 ) 
(2.7) 


The construction methods of the Riemann-Green function are described, for example, in 
[2, 13-17]. 


D. Implicit solution of the problem 

It is convenient, to write the density of a conservation law, i.e. the function ^(RRR 2 ), 
in the form (^(R 1 , R 2 ^ 1 , r 2 ) 

<p(R\R 2 \r\r 2 ) = M(r l ,r 2 )$(i? 1 , R 2 ^ 1 ,r 2 ), M(r\r 2 ) = 2 . (2.8) 

A 2 (r 1 ,r 2 ) - ARr 1 ,?- 2 ) 

The solution of (2.1), (2.2) can be represented in implicit form as [2] 

R}(x,t) = r 1 (6) = Rl(b), R 2 (x,t) = r 2 (a) = RjRa), (2.9) 
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where a, b are the new variables (Lagrangian variables). 

The connection between the new variables a, b and old variables x, t has the form 

t = t(a,b ), x = x(a,b). (2.10) 

Function t = t(a, b) is calculated using the density of the conservation law <p(R 1 , R 2 ^ 1 , r 2 ) 
and the initial data Rq(x), Rq(x) [1, 2] 

b 

t(a,b)=t 0 + ^ J ip(R 1 0 (T),R 2 0 (T)\r 1 (b),r 2 (a))dT. (2.11) 

a 

Function x = x(a, b) is calculated by analogy [2], Note, that this function is not required 
for further. We assume that this function is the given function. 

If the equations (2.10) are solvable explicitly 

a — a(x,t), b — b(x,t), (2.12) 

then we have explicit solution 

R 1 (x, t) = Rl(b(x, t)), R 2 (x,t) = Rl(a(x,t)). (2.13) 

In principle, one can assume that the original problem is solved. There is a system 
of nonlinear transcendental equations (2.10) for variables a, b , where the functions t(a,b), 
x(a, b) are completely determined. Solving this system for each fixed point (x, t ) we obtain 
the solution in the form (2.12). If we have good numerical algorithms for solving systems of 
transcendental equations and good initial approximations, then the solution of the Cauchy 
problem for ODE’s is not required. 

E. Solution on isochrones 

To construct the solution in the form (2.9) we proposed [1] to solve the Cauchy problem 
for ODE’s. We fix some value t = £*, specifying the level line (isochrone) of function t(a,b) 

t* = t(a,b). (2-14) 

We assume that the isochrone is determined on the plane (a, b) by the parametrical equations 

a = a(r), b = 6(r), (2-15) 
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where r is the parameter. 

We choose the values a*, fe* which determine some point on isochrone t — 


t* = t(a*,b*). (2.16) 

In practice, the values of a*, 5* one can choose using the line levels of function t(a,b ) for 
some ranges of parameters a, b. 

The coordinate x on isochrone, obviously, is determined by the expression 

x — x(a(r),b(r)) = X(r). (2-17) 

To determine the functions a(r), b(r), X(r) we have the Cauchy problem [1] 

= ~t b (a, b), ^ = t a (a,b ), (2.18) 

dr dr 

dx 

— = (X 2 (r 1 (b), r 2 (a)) — X 1 (r 1 (b), r 2 (a)))t a (a, b)t b (a, b), (2.19) 

a | T=0 = a*, b\ r=Q = 6*, X| r=0 = X*. (2.20) 

Here the values a*, fe* are given. To determine A"* we need to solve the problem 
dY(b) 

db =x b (a*,b) = A 2 (r 1 {b),r 2 (a*))t b (a*,b), Y(a*) = a*. (2.21) 

Integrating from a* to fe* we get 

X, = Y{K). (2.22) 

Note, that AT = x(a*, b *) is the x coordinate corresponding to r = 0. 

Integrating the Cauchy problem (2.18)-(2.20) we obtain the solution on isochrone 

E}(x, £*) = Rl(b(r)), R 2 (x,t*) = Rl(a(r)), x = X(r). (2.23) 


Moving along isochrone, that is, changing the parameter r, we obtain the solution which 
depends on x as the fixed time moment t — t*. Pay attention that the right hand sides 
of differential equations, in particular, t a (a,b), t b (a,b ) are easily computed with the help of 
(2.8), (2.9), (2.11). 

We present some auxiliary notation and relations that are useful in calculations 

v{T\a,b) = ^(r),^(r)|r 1 (6),r 2 (a)). (2.24) 


6 



Then relation (2.11) has the form 


t(a,b)=t 0 + / ip(r\a,b)dT. 


(2.25) 


The right hand sides of equation (2.18) can be easily computed 

b b 

t a (a,b) = —ip(a\a,b) + ip a (r\a,b) dr, t b (a,b) — <p(b\a,b) + / <pb(r\a, b) dr, (2.26) 


where 

(p a (r\a,b) = i^(i2j(T),i2g(r)|r 1 (6),r 2 (a))r^(a), (2.27) 

<Pb{r\a,b) = ^M- R o( r ).- R o( r )k 1 ( 6 )d 2 («)K( 6 )' ( 2 - 28 ) 

We make several important notations. First, the right hand sides of differential equations 
(2.18) can be set with accuracy to an arbitrary multiplier, which essentially override the 
parameter r. In some cases, a good choice of the parameter r allows to solve the Cauchy 
problem more effectively. Second, one should not put r — x. This replacement reduce 
the number of equations and give more natural form of the solution: R 1 (x, t*) = Rl(b(x)), 
R 2 (x,t*) = Rl(a(x)). However, this option does not allow us to construct the multi-valued 
solutions, in particular, does not allow to study the breaking solutions (see [1]). 


III. ZONAL ELECTROPHORESIS 

To demonstrate the effectiveness of the proposed method we consider the solving of the 
zonal electrophoresis equations with different initial data. The results are given in terms 
of the original equations, i.e. for values of u 1 , and in terms of the Riemann invariants R l . 
This section also demonstrates the possible simplification of the method in the case when 
the Riemann-Green function is represented in the form (1.1). 

We consider the system of equations describing the process of mass transfer under the 
action of the electric field, more precisely, the zonal electrophoresis model [8-10, 12] 

TT~) =0 > u l + =0, s = m 1 +m 2 , (3.1) 

1 + s Jx \1 + S J X 

M 1 | t=0 = ‘uS(®), u 2 \ t=0 = ul{x). (3.2) 


u. 


f/fj 2 
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Here u k are the concentrations, / i k are the effective component mobilities, s is the conduc¬ 
tivity of the mixture, ul(x), Uq(x) are given functions (initial concentration distributions). 
This system written in the Riemann invariants has the form 

R] + A 1 (i2 1 , R 2 )R l x = 0, R 2 + A 2 (R\ R 2 )R 2 x = 0, (3.3) 

A 1 (i2 1 , R 2 ) = R 1 (R 1 R 2 ), X 2 (R\ R 2 ) = R 2 {R l R 2 ). (3.4) 

The correspondence with the variables is given by the expressions 


H 0 = 1 + w 1 + w 2 , Bq = + n 2 + u 1 n z + u 2 n l , C Q — yuV 2 , 


1 I ,,2 , , 1.2 I „, 2 ,, 1 


.1 ,.2 


(3.5) 


1 1 + u 1 + u 2 R 1 + R 2 /U 1 + /i 2 + u 1 ^ 2 + u 2 /i l 


R 1 R 2 


n l n 2 


R l R 2 


fi 1 n 2 


x /i 2 (i? 1 — /j 1 )(R 2 — /i 1 ) 2 / n 1 (R 1 — /i 2 ){R 2 — n 2 ) 


R l R 2 (/i l -n 2 ) ’ R l R 2 (/i 2 -/i 1 ) 

The initial data for (3.3), (3.4) are written with the help of (3.2) as 


(3.6) 


fl'Lo = -RoM. fl2 Ln = fl oM' 


t =0 


(3.7) 


A. The Riemann—Green function and implicit solution 


In the case (3.3), (3.4) the equations (2.3), (2.4) have the form 


$ R 1 R 2 


R 1 


-$ R 1 — 


R 2 


; f h n 2 = 0. 


R 2 (R 1 -K 2 ) n R 1 (R 1 - R 2 ) 

The Riemann-Green function for the equation (3.8) is well known (see, i.g. [13-15]) 

02 , 1 2 \ _ ((R 1 + R 2 ){r l + r 2 ) — 2R 1 R 2 — 2r 1 r 2 )r 1 r 2 

1 ’ 1 ’ J R 2 R?(r'-r 2 ) 2 

It is obvious that ^(it 1 , r 2 ) can be represent in the form 


(3.8) 


(3.9) 


§{R\ R 2 \r\r 2 ) = V k (R\ R 2 )Q k (r\r 2 ), 


k =1 


where 


R 1 + R 2 


(r 1 — r 2 ) 2 ’ 


Vi(R\R 2 ) 


R l R 2 ’ 


Qi{r l ,r 2 ) 















v 2 {r\r 2 ) = 


Q 2 (r 1 ,r 2 ) = 


2 (r 1 r 


1 1 V , 2\2 


R 1 R 2 ' 

V 3 {R\R 2 ) = 1, Q 3 (r\r 2 ) = 


(r 1 — r 2 ) 2 ’ 
2r 1 r 2 


( r i - r 2 ) 2 ’ 


This representation allows to considerably simplify the calculation of the integrals in 
the expressions (2.25), (2.26). In fact, the integrals are calculated only from functions 
V k (Rl(r),Rl(r)) that contain only the initial data and do not contain the variables r 1 , 
r 2 . Moreover, for problem (3.1), (3.2) the functions Vi(R\{r), Rl(r)), V 2 (Rl(r), -Rq(t)) are 
linear combinations of the initial data Uq, Mq (see (3.5)). 

Omitting the cumbersome calculations, we present the final form of the implicit solution 
for (3.3),(3.4) 


where 


t = t(a , b) = t 0 + M w {b -a) + M tt F + MG, t 0 = 0, 
x = x(a, b) = °^ h - + M x0 (b - a) + M xF F + M xG G, 


F(a,b) = / f{r) dr, /(r) = 


Rl{r) + Rl{r) 
Rl(r)R 2 0 (r) ’ 


(3.10) 

(3.11) 

(3.12) 


G(a,b) = / g{f) dr, g(r) = 


Rl(r)R 2 o(ry 


rtF rl + r2 


M = 


t r 1 _ r 2\3 ' 


<-) 1 2 O 

M tG = , T . M t0 — 


f.pl — .^ 2^3 - 


r l _ r 2')3 J 


m , f = 2 (rV 2 ) 2 ^J^+r 2 )^ 2 ) 2 


(" r l _ r 2'\3 


. mx = -gm- 


o 3(r 2 ) 2 r 1 — (r 2 ) 3 + 3(r 1 ) 2 r 2 — (r 1 ) 3 


2^3 ’ 


M xu = 


2 (r 1 — r 2 ) 3 

R l {x,t) = r 3 (5) = Rq(P), R 2 (x,t ) = r 2 (a) = i?o( a )- 


(3.13) 

(3.14) 


(3.15) 


We recall that for solving of the problem on the isochrones it is sufficient to know only the 
function t(a,b ) which is determined by the formulae (3.10), (3.12), (3.13). The expressions 
(3.11), (3.14) are given for completeness. 
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B. The solution on isochrone 


To determine the solution R 1 (x,t Jt ), R 2 (x, i*) on isochrone £* = i(a*,6*) we have the 
Cauchy problem for the variables a, 6, F, G, A" instead of (2.18)-(2.20) 

Xr, (3.16) 

| T =0 "I — " l-r=n — " *’ ^l-r=n~ "I T=n ~ X*- (3-17) 


da 

db 

dF „ 

dG 

= G r , 

dX 

i t b , 

dr 

dr 

t — — w 

dr 

dr 

dr 

£ 

* 

cr 

II 

o 

= h, 

F Lo = F *> G \ 

T—0 

= G*, 

X\ 

1 T= 


Here, the right hand sides of ODE’s has the following form 


t a = t(a , b, F(a, b),G(a , b)), t b = j'-t(a , b, F(a, b ), G(a, b)), (3.18) 

da db 

F T (a, b) = f(b)t a + f{a)t b , G r (a, b) = g(b)t a + g(a)t b , (3.19) 

X T = (A 2 (r\b),r 2 (a)) - A \r\b),r 2 (a)))t a t b , (3.20) 

r 1 = r 1 (6) = Fq(6), r 2 = r 2 {a) = R 2 Q {a). (3.21) 


Note, that all the right hand sides of the differential equations are calculated using the 
explicit formulae (3.12)-(3.15) and the initial data (3.7). The function t a , t b do not contain 
integrals. These integrals are replaced by new variables F and G which are determined by 
the solution of the Cauchy problem. 

To determine the values F*, G*, A"* we solve the problem 

dY }!^ = (A 2 (r 1 (6), r 2 (a*))4(a*, 6, F(a*, b), G(a*, b)), (3.22) 

db 

dF(a*,b) dG(a*,6) ... 

Y(a*) = a*, F(a*,a*)=0, G(a*,a*) = 0. (3.23) 

Integrating the Cauchy problem (3.22), (3.23) from b = a* to b = 6* we get the values 
AT = Y(b *), F* = F(a*,6*), G* = G(a*,6*). 

Strictly speaking, the Cauchy problem (3.22), (3.23) just allows to calculate the integrals 
F(a*,6*), G(a*,5*) and to determine the correspondence between the parameter r and a: 
coordinate. If the explicit expression (3.11) for the function x(a,b ) is given then one can 
exclude the equations for A(r) from the Cauchy problem. 

In the next sections, the numerical solution of the problem (3.1), (3.2) or problem (3.3), 
(3.4), (3.7) with the different initial data are presented. 
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1. The periodic initial data 


We assume that the concentration of u l as time t — 0 is periodic in space 

ul(x) = u\ + 7 1 cos Q 1 x, u 2 q(x) = u 2 + y 2 sin Q 2 x, (3.24) 


where u\, y*, IT are constants. 

Physically these initial concentration distributions correspond to a periodic perturbation 
of constant concentrations ul, u 2 with amplitudes 7 1 , y 2 and periods 27r/f2 1 , 2-k /VL 2 . 

The results of calculations are given for parameters 


1 . 0 , 

p 2 = 3.0, 

ul = 1 . 0 , 

ul = 4.0, 

0.3, 

y 2 = 0.9, 

n 1 = 1 . 0 , 

Q 2 = 2 . 0 . 


(3.25) 


The initial concentration distributions and the corresponding distribution of the Riemann 
invariants are shown on Fig. 1. 




FIG. 1: Initial concentrations Uq(x), Uq{x) and the Riemann invariants R^(x), Rq(x) 

The results of calculations of the functions w 1 (a:, t), u 2 (x,t ) are shown on Fig. 2-5 as 
t = 1.744, 2.944, 3.536, 6.174, 8.843, 10.137, 14.728, 21.120. The red lines correspond to the 
initial periodic distribution, and the blue lines correspond to the concentration distribution 
at the appropriate time. For all calculations we use parameter step St = 0.05. We recall 
that the values a*, 6 * are correspond to isochrone t* = t(a*, 6 *) . 
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FIG. 2: (left) a* = —0.29, 6 * = 0.90, 5t = 0.05, x(a*, 6 *) = 1.106, £(a*, 6 *) = 1.744; (right) 
a* = —1.09, 6 * = 0.90, 5t = 0.05, x(a*, 6 *) = 1.260, f(a*, 6 *) = 2.944 



FIG. 3: (left) a* = —1.29, 6 * = 0.90, dr = 0.05, x(a*, 6 *) = 1.332, t(a*, 6 *) = 3.536; (right) 
a* = —1.89, 6 * = 0.90, dr = 0.05, x(a*, 6 *) = 1.623, f(a*,i>*) = 6.174 

In a real physical situation concentrations n 1 (x, t), u 2 (x, t ) are one-valued functions of the 
variable x. Thus, the real situation corresponds to time before the moment of the breaking 
concentration profiles (approximately, t « 3.536, the first three hgires). The remaining 
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figures illustrate the possibilities of the method which allows to construct the multi-valued 
solutions. 



FIG. 4: (left) a* = —6.89, 6 * = 0.90, 5t = 0.05, x(a*, 6 *) = 1.875, f(a*, 6 *) = 8.843; (right) 
a* = —2.89, 6 * = 0.90, 5t = 0.05, x(a*, 6 *) = 2.004, f(a*, 6 *) = 10.137 



FIG. 5: (left) a* = —7.89, 5* = 0.90, 5t = 0.05, x(a*,6*) = 2.580, f(a*,6*) = 14.728; (right) 
a* = —9.89, 6* = 0.90, dr = 0.05, x(a*,6*) = 3.228, t(a*,6*) = 21.120 

The evolution of the Riemann invariants R 1 (a:,t), R 2 (x,t ) corresponding to Fig. 2-5 is 
shown on Fig. 6-9. Note that changing of the Riemann invariants R 1 , R 2 over time is more 


13 





















‘regular’ than changing of concentrations u 1 , u 2 . 



FIG. 6: (left) a* = —0.29, 6* = 0.90, 5r = 0.05, x(a*,6*) = 1.106, f(a*,6*) = 1.744; (right) 
a* = —1.09, 6* = 0.90, 5t = 0.05, x(a*,6*) = 1.260, i(a*,£>*) = 2.944 



FIG. 7: (left) a* = —1.29, 6* = 0.90, 5t = 0.05, x(a*,6*) = 1.332, t(a*,6*) = 3.536; (right) 
a* = —1.89, 6* = 0.90, 5t = 0.05, x(a*,6*) = 1.623, i(a*,£>*) = 6.174 
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FIG. 8: (left) a* = —6.89, 6* = 0.90, 5t = 0.05, x(a*,6*) = 1.875, t(a*,6*) = 8.843; (right) 
a* = —2.89, 6* = 0.90, <5r = 0.05, x(a*,6*) = 2.004, t(a*,6*) = 10.137 



a : 


9.89, 6* = 0.90, 5t = 0.05, x(a*,6*) = 3.228, f(a*,6*) = 21.120 














2. Wave packet 


We assume that the concentrations u l as time t — 0 are 

ul = ul(l + ^e~P lx2 cosQ}x), (3.26) 

Up = w 2 ( 1 + 7 2 e _ ^ 23;2 sinf 2 2 x), 
where u\, 7 *, fP, f3 l are the constants. 

Physically this initial concentration distribution corresponds to the wave packet pertur¬ 
bation. 

The results of calculations are given for parameters 

pi 1 = 1.0, p 2 = 3.0, u\ = 1.0, ul = 4.0, (3 1 = f3 2 = 0.1, (3.27) 

7 1 = 0.1, 7 2 = 0.3, fl 1 = 5.0, n 2 = 5.0. 

We restrict only an illustration of the calculations for time t = 1.655, 3.737. On Fig. 10, 11 
the distribution of the concentrations and the Riemann invariants, respectively, are shown. 
At f ~ 1.655 the braking concentration profiles (and, of course, the Riemann invariants) are 

well visible. 



FIG. 10: (left) a* = —0.29, 6* = 0.90, St = 0.05, x(a*,6*) = 1.238, i(a*,6*) = 1.655; (right) 
a* = —1.29, 6* = 0.90, 5t = 0.05, x(a*,6*) = 1.673, t(a*,6*) = 3.737 
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FIG. 11: (left) a* = —0.29, 6* = 0.90, Sr = 0.05, x(a*,6*) = 1.238, t(a*,6*) = 1.655; (right) 
a* = —1.29, 6* = 0.90, St = 0.05, x(a*,6*) = 1.673, f(a*,6*) = 3.737 


5. The Gaussian distribution of the initial concentrations 


The initial Gaussian distribution of the concentrations is most often realized in zonal 
electrophoresis experiments. 

We take the following initial data 


u 


; = ul + ye-" w>", ul = ul + 


(3.28) 


where u\, 7 *, x®, /3® are the constants. 

The results of calculations are given for parameters 

/j 1 = 1.0, /r 2 = 3.0, ui = 1.0, ul = 4.0, /3 l = /3 2 = 0.3, (3.29) 

7 1 = 2.0, 7 2 = 3.0, x 1 = 5.0, x 2 = 0.0. 

On Fig. 12 the isolines of the functions t(a, b ), x(a, b ) on the plane (a, b) are shown. 
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FIG. 12: The isolines of the function t(a,b ) (left) and x(a,b ) (right) 


On Fig. 13-16 for time moments t = 1,146, 2.338,3.438, 5.614, 6.358, 10.192, 16.350, 
22.513 the results of calculations of the functions u k {x , t) are shown. The red lines correspond 
to the initial distribution, and the blue lines correspond to the concentration distribution in 
the appropriate time. We use parameter step 5t = 0.05 for all calculations. We recall that 
the value a*, fe* are correspond to isochrone £* = i(a*, 6 *) . 



FIG. 13: (left) a* = 0.49, 6* = 0.90, 5t = 0.05, x(a*,6*) = 1.009, t(a*,6*) = 1.146; (right) 
a* = 0.09, 6* = 0.90, dr = 0.05, x(a*, 6*) = 1.122, t(a*, 6*) = 2.338 
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FIG. 14: (left) a* = —0.29, 6* = 0.90, dr = 0.05, x(a*,6*) = 1.226, t(a*,6*) = 3.438; (right) 
a* = —1.29, 6* = 0.90, dr = 0.05, x(a*,6*) = 1.435, i(a*,6*) = 5.614 



FIG. 15: (left) a* = —1.89, 6* = 0.90, dr = 0.05, x(a*,d*) = 1.509, t(a*,6*) = 6.358; (right) 
a* = —4.89, 6* = 0.90, dr = 0.05, x(a*,6*) = 1.921, t(a*,6*) = 10.192 


19 





































i—I—■—■—■—I——1—■—■—I—I—I—■—■—■—I—I—■—I—■—I—I—■—I—■—I 

-5 20 


9 ■ 




0.05, x(a*,6*) = 2.591, t(a*,b*) = 16.350; (right) 


FIG. 16: (left) a* = —8.89, 6* = 0.90, St 
a* = —12.89, 6* = 0.90, St = 0.05, x(a*, &*) = 3.262, f(a*, &*) = 22.513 


We again recall that in a real physical situation the concentrations ■u 1 (a:,t), u 2 (x,t ) are 
the one-valued functions of the spatial variable x. The real situation corresponds to time 
before breaking concentration profiles (approximately, t ~ 6.358, the first five figures). The 
remaining figures illustrate the possibilities of the method which allows to construct the 
multi-valued solutions. 

The Fig. 17 shows the results of calculations for the case when the initial data have the 
form 

ul = ul, ul = ul + ye - ^ 2 , (3.30) 

H 1 = l.o, /i 2 = 3.0, ul = 0.5, ul = 4.0, 0 = 0.3, 7 = ±1.0. 

The results illustrate that the presence of the initial perturbation of any only one concen¬ 
tration (in this case u 2 ) leads to the perturbation of the other concentrations (in this case 
u l ). Strictly speaking, the result is quite obvious, as it is easy to check that the original 
equation (3.1) does not admit the solutions u 1 = const, u 2 7 ^ const. 

It is interesting to note that the initial Gaussian’s perturbation single peak, in the process 
of evolution, splits into two peaks. For initial data (3.30), if 7 = +1, the distribution of 
perturbations occurs quite regularly, whereas in the case 7 = — 1 , over time there are the 
breaking concentration profiles, and the occurrence of shock waves. 
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FIG. 17: (left) 7 = +1, a* = -8.89, b* = 0.90, 5t = 0.05, x(a*,6*) = 4.268, i(a*,6*) = 18.091; 
(right) 7 = —1, a* = —8.89, 6 * = 0.90, 5t = 0.05, x(a*,6*) = 6.467, t(a*,6*) = 19.942 


IV. CONCLUSIONS 

Implicit solution (3.10), (3.15) of the original problem (2.1), (2.2), of course, is very 
important in studying of the original Cauchy problem properties. However, from our point 
of view, the implicit solution is not less complex than the original problem. For practical 
applications we must get, in one way or another, the explicit relations (2.12). In the general 
case we need to use a numerical methods, for example, Newton’s method, for solving of 
the transcendental equations. This, in turn, requires a good initial approximations, or the 
using of the movement parameter method. It is especially difficult to numerically solve 
the system (2.10) in the case when original problem has multi-valued solition. In other 
words, using of the numerical methods of solving the systems of transcendental equations 
is not much easier than the application of the direct methods for solving of the original 
problem, for example, using finite difference or finite volume method. Instead solving of 
the transcendental equations system we propose to solve the Cauchy problem for ODE’s. 
Even if it is need to solve ODE’s numerically, for example, by the methods of Runge-Kutta 
method, the numerical algorithm is realized simpler than the algorithm for solving system 
of nonlinear transcendental equations. 
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